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LHCb, bolstered up by the 10 12 Z?-hadrons produced yearly, is an excellent place to study rare 
B decays. The B®—t fJ. + jJ,^ decay, generated by flavour-changing neutral current, is strongly 
suppressed within the Standard Model (SM). However, its branching ratio can be significantly 
enhanced if New Physics (NP) exists. The current best limit is still an order of magnitude above 
the SM prediction, thus leaving room for observation of NP effects. In this document, the current 
LHCb analysis strategy for the B® — ► search is presented. 
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1. Introduction 

The Standard Model (SM) is successful in explaining almost all observations in particle physics 
experiments so far. Nevertheless, there are reasons to consider it as a low energy effective limit of 
a more general theory In that prospect, observables for processes where the SM contribution is 
highly suppressed are particularly interesting. Within the SM, flavour-changing neutral current 
processes are highly suppressed since they are forbidden at tree level and can only proceed via 
loops diagrams. If New Physics (NP) exists, new particles can contribute to those processes in the 
loop diagrams, modifying observable quantities with respect to the SM prediction. The B®— > jU + /i~ 
decay branching fraction is such an observable. LHCb ^ will take advantage of the copious b- 
hadron expected at LHC [gj. LHCb is a single-arm forward spectrometer primarily optimized to 
the study of CP- violation and rare decays in ^-hadrons. The detector is characterised by its precise 
vertex detector, powerful particle identification capabilities and versatile trigger. Nominally, LHCb 
will operate at a luminosity ££ = 2 x 10 32 cm~ 2 s _1 , giving 2 fb _1 per year (10 7 seconds) of data. 
The analyses presented in this document are applied to Monte Carlo simulated data with a full de- 
tector response, including pile-up (multiple pp collisions in a single bunch-crossing) and spill-over 
(signal coming from particles produced in a previous bunch-crossing). 

Within the SM, the B® —> /J. + H~ decay occurs through loops diagrams like the one in Figure 
and its branching ratio is expected to be BR(B° -> jU+/j. ) = (3.35 ± 0.32) x 1(T 9 @|. In minimal 
supersymmetric extensions of the SM (MSSM), this decay would receive additional contributions 
from diagrams of the kind shown in Figure f| The branching ratio is then proportional to the 
sixth power of the ratio of the Higgs vacuum expectation values, tan j8 , and can be considerably 
enhanced [Q]. The current upper limit given by the two Tevatron experiments is 4.7 x 10 at 90% 
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Figure 1: Feynman diagram of the Z?' 
decay within the SM. 
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Figure 2: Feynman diagram for a MSSM con- 
tribution to M + jU~ . 
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2. Strategy for the 5" — > Decay Search 

The current LHCb strategy for the — > £i + jU~ decay search [j|] can be summarized as fol- 
lows: 

• The analysis relies on three independent variables, related to the dimuon invariant mass, the 
particle identification of the daughters, and geometrical information from the decay topology, 
respectively. Since the three variables are uncorrected, they can be calibrated independently. 
The calibration methods are designed to rely solely on real data. 

• For each event, likelihood values are calculated for the each of three variables, under the 
signal and background hypothesis. 

• The compatibility of the obtained distributions of the likelihood values is tested against 
B® — > H + IJ.~ branching ratio hypotheses, using the CL modified frequentist method [^] with 
the calibrated likelihood distributions for signal and background. The final result is either a 
measurement or an upper limit of the branching fraction. 

• Since the number of 5° produced is not precisely known, the use of a normalisation channel 
with a well known branching fraction is required to obtain an absolute measurement of the 
branching fraction or upper limit. 

2.1 Invariant Mass Likelihood 

The signal likelihood distribution is calibrated using the B® — > K + K decay, which is kine- 
matically very close to the signal. The selection of the B® — > K + K events uses a cut on the kaon 
identification, which biases the B° s mass resolution. The effect of the kaon selection cut on the 
mass resolution is assessed onB?>-> h + h~ control channels, and the obtained correction factor is 
applied to B® — ► K + K sample (see Figure ||). 

The background likelihood is calibrated using B° s — ► ;U + /i~ candidates in the mass sidebands, 
outside the ±60MeV/c 2 region around the nominal B® mass. 

2.2 Particle Identification Likelihood 

The identification of the daughter particles as muons relies mainly on the muon system in- 
formation. The signal likelihood is calibrated using a high-purity muon sample obtained from 
J/\l/—>- M + jti _ decays where only one muon is identified. 

The background likelihood due to misidentified hadrons is obtained from A— >p %~ decays, 
which are very abundant in LHCb and have a very clean signature. 
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M (MeV/c 2 ) 

Figure 3: Invariant mass distributions for B° s ^> K + K~ with a correction for the kaon identification cut 
effect (solid curve) and without correction (dashed curve). Points with error bars represent B® — ► jl + ix~ 
invariant mass distribution. 



2.3 Geometrical Likelihood 

The geometrical likelihood is given by a single real number computed from the following 
topological variables: 

• the B Q S impact parameter with respect to the primary vertex, 

• the proper time, 

• the smallest impact parameter significance of the two muons with respect to any primary 
vertex, 

• the distance of closest approach between the two muons, 

• the isolation of each muon track. 



Using Monte Carlo samples, this likelihood is defined to take values between and 1, and to 
have a flat distribution for signal events and a distribution peaked at for background events, as 
shown in Figures |] and ||. 

The signal distribution will be calibrated by real data using ZJ9 — > h + hT events, which are 
topologically identical to B® — > /i + jU~ and the background distribution by events in the B® mass 
sidebands. 



3. Normalisation 

The branching fraction is expressed as: 



2 Ofofo =S?jnt X fs X £ s ig 
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Figure 4: Geometrical likelihood distribution of the 
signal obtained by simulation. 




Figure 5: Geometrical likelihood distribution of 
background events containing two muons in the de- 
tector acceptance. 



where <y^i g is the number of signal events, o bb - is the bb production cross-section, J2f mt the integrated 
luminosity, f s the probability of a b quark to hadronize in a B° s meson, and £ s j g the product of the 
reconstruction, trigger, and selection efficiencies. Since c bb - JZi nt will not be precisely known, 



we measure the B; 



o 



jU + jU branching fraction using a normalization channel whose branching 



fraction is well known. The B® — ► jU + jU branching fraction is then: 



BR(fl°->/i + /i _ ) = BR n 
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where / nor m> £norm and jV nwm are the quantities for the normalization channel, with definitions 
analogous to those of the signal channel. 

Possible normalisation channels are B + — ► J/y(I-LI^)K + and B° — ► K + n~ . Particular care has 
been taken to analyze the signal and normalization channels in a common way so that any large 
systematic effects in the efficiency ratio cancel. 

The main systematic arises from the ~ 13% uncertainty on the ratio -f- or With sufficient 

Is Js 

statistics, all other systematic uncertainties are expected to get much smaller than that, as the anal- 
ysis relies solely on data. In the future, pending on a more precise measurement of its branching 
fraction at Belle, B® — > D~ n + can be used as normalization channel, with the effect of eliminating 
the use of the factor 

Is 



4. Analysis Sensitivity 

The expected sensitivity of LHCb to the B® — > jU + jU~ decay as a function of the integrated 
luminosity is shown in Figure |6| The solid line in the left plot shows the expected upper limit, 
at 90% confidence level, on the branching fraction when no signal is observed, for pp collisions 
at = 8 TeV. The branching fraction for which a 5a discovery or for which a 3a evidence is 
expected, is shown on the right plot for 14 TeV collisions. 

The left plot also indicates the expected final limit from the Tevatron experiments, extrapolat- 
ing the current results to 8 fb _1 of data per experiment. It shows that LHCb competes with the 
current Tevatron limit (4.7 x 10~ 8 ) with approximately 0.1 fb~ ! of data and overtakes the final 
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Figure 6: Expected 90% CL upper limit on BR(Z? S — * ) in absence of signal for 8 TeV collisions (left) 

and BR(Bj — > ) at which a 5(7 discovery, stars, or a 3(7 evidence is expected, plain blue curve, for 14 

TeV collisions (right) as a function of the integrated luminosity. On the exclusion plot, the black curve is the 
result of the standard analysis, and the red dashed curve is the result of the robust analysis. The background 
estimate has conservatively been set to its 90% CL upper limit in the exclusion plot, and similarly the dashed 
curves indicates the 90% CL upper and lower limit in the observation case. The Tevatron limit is calculated 
by extrapolating the current result to 8 fb _1 per experiment. 

expected limit with about 0.2 fb _1 . NP models with high tanf} value are strongly constrained in 
the process. In that time frame, the detector may not be fully understood yet. Therefore, an alter- 
native robust analysis is designed using variables with a similar physical content to the ones of the 
standard analysis, but avoiding the use of error estimates. This implies a modified selection and 
definition of the geometrical likelihood. The robust analysis sensitivity is also depicted in Figure || 
(left), as a red dashed curve. The robust analysis presents a sensitivity compatible with the one of 
the standard analysis. Therefore, it constitutes a valuable option for the early data. 

About 3 fb^ 1 are enough for a 3a evidence if the branching fraction is the SM prediction. 
Any enhancement driven by NP will be observed sooner. Particularly, if the branching fraction 
is as high as 2 x 10~ 8 , as predicted in Reference [g], a 5a discovery is possible with very little 
luminosity (< 0.4 fb _1 ). With 10 fb _1 , a 5a discovery occurs if the branching fraction is at the 
level of the SM prediction. 
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